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Thermal fluctuations can play an important role in the buckling of elastic objects at small scales, 
such as polymers or nanotubes. In this paper, we study the finite-temperature buckling transition 
of an extensible rod by analyzing fluctuation corrections to the elasticity of the rod. We find that, 
in both two and three dimensions, thermal fluctuations delay the buckling transition, and near the 
transition, there is a critical regime in which fluctuations are prominent and make a contribution 
to the effective force that is of order ^/T. We verify our theoretical prediction of the phase diagram 
with Monte Carlo simulations. 
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I. INTRODUCTION 

When a thin elastic rod is under compression on its two 
ends, it experiences an instability towards buckling as the 
compression exceeds a critical value; this is the classical 
Euler buckling problem [1, 2]. This critical compression 
is determined by the competition between the compres¬ 
sion and bending energy costs of the rod. The buckling 
instability plays an important role in many problems in 
fields ranging from physics to engineering and biology [3- 
8 ]. 

More recently, experimental studies on buckling phe¬ 
nomena at small length scales, such as the buckling of stiff 
or semiflexible polymers, nano-filaments, and nanotubes, 
have been enabled by advances in various technologies [9- 
13]. These studies may lead to novel devices that utilize 
transitions between multiple mechanical ground states. 
At these small scales, it is necessary to include effects 
of thermal fluctuations, which have been shown to lead 
to interesting phenomena near mechanical instabilities 
in various systems [14-20]. Such thermal-fluctuation 
effects have been theoretically investigated, and phe¬ 
nomena such as corrections to the projected end-to-end 
length, shifts in the critical compression, and softening 
of the buckling transition have been discovered [21-28]. 
However, most of these theoretical studies on how fluctu¬ 
ations renormalize the buckling transition have focused 
on the case of inextensible polymers and have employed 
the worm-like chain model, which assumes that the poly¬ 
mer has a constant contour length. This is an idealized 
limit where the rod cannot be stretched/compressed. For 
real rods, although the resistance against stretching is 
much stronger than that against bending, it is worthwhile 
to discuss whether the extensibility of the rod changes 
what is known about buckling at finite temperature. 

In this paper, we investigate finite-temperature buck¬ 
ling using a model elastic energy that allows for rod ex¬ 
tensions. In this model, the end-to-end distance is the 
control parameter (fixed-strain ensemble), and the rod 
is allowed to have transverse fluctuations, which both 
stretch/compress and bend the rod. By integrating out 
higher-momentum modes which couple to the first fun¬ 


damental mode through anharmonic terms, we calculate 
fluctuation corrections to the rigidity and analyze the 
buckling transition of the renormalized theory. We find 
that, in both two and three dimensions, thermal fluctu¬ 
ations shift the buckling transition to larger-magnitude 
values of compression. Our Monte Carlo simulations ver¬ 
ify the analytic phase diagram we obtain (Fig. 1). In 
addition, we also analytically calculate the effective force 
of the rod, showing that, close to the buckling transi¬ 
tion, thermal fluctuations are prominent and contribute 
an 0{VT) correction to the effective force. 

It is worth pointing out that, in the presence of thermal 
fluctuations, the rod is never completely “straight.” The 
physical meaning of having a “straight-buckled” transi¬ 
tion is that the mean-square transverse fluctuations of 
the rod (e.g., the mean-square transverse displacement of 
the midpoint of the rod) change from zero in the straight 
phase to a nonzero value in the buckled phase. In other 
words, the elastic free energy minimum of the rod changes 
from the straight configuration to the buckled configura¬ 
tions. 

This paper is organized as follows: we construct the 
model and discuss the analytic theory in Sec. II and 
present the Monte Carlo simulations in Sec. III. Then, in 
Sec. IV, we summarize our results and discuss relations 
to other studies. 


II. MODEL AND ANALYTIC THEORY 

A. The extensible-rod Hamiltonian with 
anharmonic terms 

We consider a thin elastic rod with rest length Lrest 
embedded in d dimensions, as shown in Fig. la. Here, 
Lo is the end-to-end distance, or projected length, of the 
rod, which is the control parameter of our theory, and 
L is the instantaneous contour length in the presence of 
thermal fluctuations. 

Assuming that the rod is made of a homogeneous ma¬ 
terial with Young’s modulus E, its stretching rigidity g 
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FIG. 1. (Left) Illustration of an extensible rod under compression. (Right) Predicted phase diagram for two-dimensional 
finite-temperature buckling of an extensible rod in the plane of normalized compression and temperature (as defined in Sec. II). 
The thick black curve represents the phase boundary between the straight (unshaded) and the buckled (shaded with vertical 
lines) phases. The dashed curve denotes the boundary of the critical regime (light blue region) where the thermal-fluctuation 
correction to the effective force is of 0{VT). 


and bending rigidity n are given by 
g = Ira Ejl/Yesti 

K = El = {n/A)Ea‘^, (2.1) 


where Xd parametrizes the rod using the projected end- 
to-end distance, r{xd) is the position of the rod at Xd, 
and xj_ denotes the transverse displacement of the rod. 
We define derivatives 


where a is the radius of the rod, and I is the mo¬ 
ment of inertia of the cross-section. Throughout this 
paper, we require that the relative strengths of the (me¬ 
chanical) rigidities against bending and against stretch¬ 
ing/compression of the rod satisfy 


k± K 



( 2 . 2 ) 


meaning that it is much more energetically costly to 
stretch/compress the rod than it is to bend the rod. This 
is satisfied by most microscopic rod-like objects, includ¬ 
ing polymers, nanowires and nanotubes [27]. 

The instantaneous stretching/compression elastic en¬ 
ergy of the rod can be written as 

Use — Tj-est) 

= ^g[(L-Lo) + {Lo-L,,st)]^. (2.3) 

We define r to be the force applied to the ends of the 
straight rod at T = 0, when there are no thermal fluctu¬ 
ations (i.e., L = Lq): 


— 9{Uo .^rest)5 (^■'^) 


so that T > 0 corresponds to stretching of the rod, while 
T < 0 corresponds to compression. 

We proceed to derive the Hamiltonian of the rod for 
a given compression r and an instantaneous fluctuation 
configuration, which is described by 


r{xd) = {xi,..., Xd-i,Xd) = {x±,Xd), 

0 <Xd < Lq, (2.5) 


n{xd) = 


C?X_L 

dxd 




( 2 . 6 ) 


where tt is a one- (two-) dimensional vector for the case 
of a rod embedded in two (three) dimensions. 

The Hamiltonian can then be written as an expansion 
up to 0(7r'*), 


H^Ho + H2 + H4, (2.7) 


where 

Ho = ^ (2.8) 

is the energy of the straight rod with no fluctuations [L = 
Lq): 


H2=^ 


dxd {t\tv\^ + kI'^T) 


(2.9) 


contains terms quadratic in tt, and 


1 


rLo 


= 2 y dxd 

Lq 


- ykl^lTTk - ■ tt')^ 


JJ dxddx'a\Tv{xd)\‘^\TT{x'^)\‘^ ( 2 . 10 ) 


includes terms quartic in tt. Here tt' is shorthand no¬ 
tation for d'K/dxd- This Hamiltonian iJ includes con¬ 
tributions from both stretching/compression as well as 
bending of the rod, and the details of its derivation are 
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included in App. A. The last term in H 4 , coming from 
g{L — /2, appear to be nonlocal; however, as we shall 

see, it simply leads to a term in Fourier space with 
its momentum sum limited to a special channel. 

Note that this formulation with fixed end-to-end dis¬ 
tance is the same as the one used in the classical Euler 
buckling problem in textbooks [2]. A similar formulation 
has also been used in Refs. [23, 24], which focus on quan¬ 
tum aspects of buckling. Additionally, because we are 
interested in the case where the rod is much more resis¬ 
tant to stretching than it is to bending, the stretching can 
be taken to be small and highly homogeneous throughout 
the rod. This allows for the approximation to be made 
that the parameters k and g are uniform along the rod, 
as in Ref. [29]. 


B. Classical (T = 0) Euler buckling 

The T = 0 buckling transition is obtained by analyzing 
the stability of the quadratic coefficient of the Hamilto¬ 
nian H, while the T = 0 configuration is determined by 
the location of the minimum of H. It is convenient to an¬ 
alyze this Hamiltonian in momentum space. In order to 
use the convenient exponential form of the Fourier trans¬ 
form, we employ the trick of extending the end-to-end 
distance of the rod to Xd € [—Lq,Lo] to obtain periodic 
boundary conditions from the physical fixed boundary 
conditions that x_l = 0 at = 0, Lq (further discussion 
of this can be found in App. A). The quadratic-order 
Hamiltonian, which is sufficient to ascertain the stability 
of the system, can then be written as 

8L0 Y 

where 

g=^, neZ\{0}. (2.12) 

bo 

Although the sum seemingly counts excess modes by in¬ 
cluding both positive and negative values of g, these 
modes are not actually independent: because Tr(xd) is 
real and even, satisfies the constraints that 

TTq=TT_q=TT*, (2.13) 

so that the above sum is even in q, and the number of 
independent modes is the same as in the case of expand¬ 
ing H in terms of sm{mTXd/L q). It is straightforward 
to extract the T = 0 Euler buckling condition from this 
equation. The magnitude of the lowest allowed momen¬ 
tum mode is qi = tt/Lq, since the q = 0 mode is excluded 
by the above fixed-end boundary conditions. In order for 
the Hamiltonian to have a stable equilibrium at = 0, 
its matrix representation must be positive definite - all 
its eigenvalues must be positive: 


Applying this condition to the lowest mode, we obtain 
the critical compression 


Tc(0) = (2.15) 

bo 

where the 0 in parentheses indicates that this is a T = 0 
result. Recall that t < 0 corresponds to compression of 
the rod, so that for any compression t > Tc(0) (i.e., com¬ 
pression with a magnitude less than that of the critical 
value), the rod remains straight. 

For T < rc(0), on the other hand, the harmonic-level 
Hamiltonian is no longer stable at tt = 0. The number of 
modes that have become unstable depends on the value 
of r; for (n + I)^Tc(O) < t < n^Tc(O), the first n modes 
are unstable, as each of their coefficients in H 2 is neg¬ 
ative. Thus, for increasingly negative values of r, it is 
possible to have various metastable states corresponding 
to higher orders of buckling; for any value of t < rc(0), 
however, the most energetically favorable buckled config¬ 
uration is the n = 1 mode. In this paper, we will only 
be concerned with analyzing the instability of the first 
momentum mode when considering the buckling transi¬ 
tion; therefore, our discussion in the buckled phase will 
be restricted to the case of 4tc( 0) < r < Tc{0), since 
the second mode becomes unstable for t < 4rc(0). In 
this range of compression values, the new stable state - 
corresponding to the n = 1 buckled phase - is fixed by 
the anharmonic terms in H 4 with only tti nonzero. As 
detailed in App. B, the stable configuration is described 

by 


|x_L(a;d)| 


f^ / 4Lo(Te(0) ^r) 

Ti" V 9 



(2.16) 


in the T = 0 buckled phase, where we have applied 
the limit of stretching stiffness being much greater than 
bending stiffness [Eq. (2.2)] to obtain gL^ |tc( 0)|, |r| 
to simplify the expression. 


C. Fluctuation corrections to stability and 
flnite-temperature bnckling 

The finite-temperature phases are determined by the 
minima of the free energy of the rod, which includes en- 
tropic contributions. At finite temperature, thermal fluc¬ 
tuations excite all modes of the rod, and these fluctua¬ 
tions renormalize the stability of the rod against buck¬ 
ling. In order to analyze this entropic effect on the buck¬ 
ling transition, at which the first mode qi becomes unsta¬ 
ble, we follow a procedure similar to that of momentum 
shell renormalization. We first separate the first modes 
from the higher-momentum fluctuation modes: 


T -f > 0 Vg. 


(2.14) 


(2.17) 
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where 




q if kl = 9i 
0 if \q\ > qi 


and 


0 


if |g| = qi 
if |g| > qi 


(2.18) 


It follows that Tv{xd) = Tt’^{xd) + TT^{xd)- These two 
components are decoupled in the quadratic Hamiltonian 
i? 2 ) but H 4 has cross terms. 

The partition function can then be written as 


Z = 


X)t^<T>TV^ g-y (ffo+^f2(7r<)+ff2(7r>)+ff4(7r<,7r>)) 


= / Vtt< e-^ 


')/T 


(2.19) 


where, in the second line, we define the Landau free en¬ 
ergy 

F<{7Z<) = Ho+H2{7T<) 

-TlnJ X> 7 r> e"T(^2(7r>)+ff4(7r<.7r>)) 

( 2 . 20 ) 

This Landau free energy, with all other modes tt^ in¬ 
tegrated out, determines the finite-temperature stability 
of the first mode, as discussed in detail in App. C. The 
resulting includes terms quadratic order in tt^, 


to the T = 0 buckling transition, we can expand A{t) in 
powers of f — 1, 

^+(S ■ 

The magnitude of the effective compression, f, as well 
as the effective bending rigidity, k, both decrease with 
increasing temperature. It is easy to understand the de¬ 
crease in f: thermal fluctuations tend to increase the in¬ 
stantaneous arc length of the rod from its T = 0 straight- 
rod length so that the rod effectively feels less compres¬ 
sion. In previous works [25, 30], the fluctuation correc¬ 
tion to k was shown to have a prefactor of {d — 2) in¬ 
stead of {id — 1) as we have here. The difference arises 
from the fact that the rod is assumed to be inextensible 
and, therefore, is modeled as a worm-like chain in these 
previous papers, whereas it is extensible in our model. 
Consequently, it was necessary to reparametrize the rod 
in terms of Xd rather than the arc length, s, modifying 
the form of the bending energy. 

The buckling transition occurs when the first mode 
becomes unstable, which is when 

T + kqf = 0. (2.26) 

This condition can be solved to obtain a critical tem¬ 
perature separating the straight (T > T),) and buckled 
(T < Tc) phases of the rod for a given compression f > 1, 

- _ T — 1 

^ [(d - l)g + {d+ l)f - (3d - 1)] Aif) • 


= + ( 2 . 21 ) 

with the original elastic parameters replaced by renor¬ 
malized ones. The renormalized parameters are given by 

f = T + A(f)'f |tc( 0 )| [{d - l)g -I- (d -I- l)f] 

k = K-A{f)f{M-l)K, (2.22) 


The phase boundary in two dimensions determined by 
this equation is plotted as the solid black line in Fig. 1. 
The three-dimensional version is shown in Fig. 2. In the 
limit that we have been considering of stretching stiffness 
much greater than bending stiffness [Eq.(2.2)], we have 
5 3> 1 , so that we can write a simplified expression for 
the critical temperature, 


fc{f) 


T — 1 4 (r — I) 

{d- l)gA{f) 3{d- l)g' 


(2.28) 


where we have defined the dimensionless quantities 
T = 

9 = 
f = 

and 


The expression following the arrow is the limiting case 
true for sufficiently low temperatures such that gT 1, 
since, as we can see from the initial equality in Eq. (2.28), 
that condition necessitates that f — I <C I, as well. In 
that case, we can write the critical temperature to leading 
order in f — 1 , allowing us to use the zeroth-order term 
in the expansion of A{f) in Eq. (2.25). 

This leading-order relation can be inverted to obtain 
an expression for the critical compression for buckling at 
a finite temperature T, 


Lo|re(C 

9 L 0 

|rc( 0 )| 


=( 0 ) 


> 0 , 


(2.23) 


Af) = Y. 

n=2 



3 r — 1 TT rz 

—7--T-;= cot TTVT 

2 (f-I)f 2 Vf 


(2.24) 

In accordance with the range of r we are considering in 
this paper, 0 < f < 4. As will be justified shortly, close 


a(T) 


a(o) 



3(d- I) 
4 



(2.29) 


This clearly represents a critical compression that is 
of larger magnitude than the zero-temperature critical 
value. In other words, the buckling transition is “de¬ 
layed” by thermal fluctuations. 
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FIG. 2. Predicted phase diagram for three-dimensional finite- 
temperature buckling of an extensible rod in the plane of nor¬ 
malized compression and temperature. The thick black curve 
represents the phase boundary between the straight (un¬ 
shaded) and the buckled (shaded with vertical lines) phases. 
The dashed curve denotes the boundary of the critical regime 
(light blue region) where the thermal-fluctuation correction 
to the effective force is of 0{VT). 


where 

Z> = J (2.33) 

is the quadratic-order partition function of tt^. The co¬ 
efficients 62 and 64 , and the integral over tti, are de¬ 
rived in App. D. We have only needed to retain terms to 
quadratic order in in (2.32) because modes are 
stable at = 0 ; quartic-order terms in (renormalized) 
TTi are necessary, however, because the quadratic-order 
coefficient, 62 , can become negative for tti - thus, higher- 
order terms in the potential are needed to evaluate the 
free energy. 

As detailed in App. D, we find that thermal fluctua¬ 
tions reduce the compressional force in the straight phase 
but enhance it in the buckled phase; these modifications 
are of order T except very close to the transition for small 
values of gT, where there is a decrease in the compression 
of order ^/T, as shown in Fig. 3. 


D. Effective force 


III. MONTE CARLO SIMULATIONS 


In this paper, we have utilized the ensemble with fixed 
end-to-end distance Lq. At T = 0, taking the derivative 
of the Hamiltonian with respect to Lq yields that the 
force on the rod is simply r = g{LQ — Lj-est) in the straight 
phase (with r < 0 corresponding to compressional force) 
and Tc(0) in the buckled phase (App. B). 

At finite T, we determine the effective force / through 


dF _ dr dF _ dF 
BLq SLq dr ^ dr' 

with the free energy given by 


(2.30) 


F = -rinZ = -Tln J 


(2.31) 


where, as defined earlier, F^ is the Landau free energy 
with only integrated out. It is useful to note that / is 
calculating by taking the derivative of F with respect to 
the compression r, rather than by taking the derivative 
directly with respect to Lq. This is intentional, as the 
derivative with respect to Lq would also act on the pref¬ 
actors of Lq in the Fourier transform (or, equivalently, on 
the integration limits in real space), which would intro¬ 
duce an ultraviolet divergence that scales linearly with 
the high-momentum cutoff. Strictly speaking, the effec¬ 
tive force / obtained via differentiation with respect to t 
describes the change of the free energy that occurs with 
changing the amount of compression r while keeping Lq 
constant. 

The Landau free energy F^', as defined in Eq. (2.20), 
can be written to leading order in T as 

F< = Hq- Tin Z>+b 2 T\TVi\^+ biT\ni\\ (2.32) 


We perform Monte Carlo (MC) simulations in two and 
three dimensions to corroborate our analytical results. 
The rod is discretized into N segments with fixed verti¬ 
cal length io = Lq/N along the Xd-axis. The segments 
are allowed to have transverse fluctuations x_lj- and to, 
consequently, cause stretching/compression and bending 
of the rod, as discussed in Sec. IT The fixed boundary 
conditions necessitate that xj_.o = = 0 - 

The Metropolis algorithm is used in our Monte Carlo 
simulations, in which, at each MC step, a segment is 
selected at random, and a random trial displacement in 
the transverse direction is attempted. For a given rod un¬ 
der a certain compression, runs are performed at various 
temperatures. We choose the transverse displacement of 
the middle segment, |x_|_ wj, to be our order parameter. 
In the straight phase |xj_ w | is governed by a Gaussian 
distribution with its mean at 0 , whereas in the buckled 
phase, the distribution of |x_|_ w | becomes double-well (in 
d = 2) or Mexican-hat (d = 3) with minima at 


1 /4Lo(t-c( 0) - r) 


(3.1) 


At the buckling transition, the distribution sharply devi¬ 
ates from Gaussian. To capture this transition, we cal¬ 
culate the Binder cumulant of the distribution [31], 


(|xx,fn 

3(|x^.«IT‘ 


(3.2) 


The value oIUl decreases as the temperature is lowered 
and the system experiences the straight-to-buckled phase 
transition. This decrease becomes increasingly sharp for 
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FIG. 3. Plots of the dimensionless effective force in both two (left) and three (right) dimensions at various values of T. Note 
that — / = //|'rc(0)|. Negative values correspond to a compressive force, while positive values correspond to a stretching 
force. The dots on each curve indicate the transition point between the straight and buckled phases for each value of T. The 
T = 0.001 curve in three dimensions does not have a dot as there is no phase transition at that temperature; the rod remains 
in the straight phase. 


progressively larger systems, and the simultaneous cross¬ 
ing of Binder cumulant curves for various system sizes 
determines the location of the critical temperature Tc- 


To verify our phase diagrams in Figs. 1 and 2 via the 
crossing of the Binder cumulant curves, we simulate rods 
containing 10, 12, and 14 segments (corresponding to 
Lq = 1.0,1.2,1.4, respectively, so that io — 0.1 is kept 
fixed). As discussed in Sec. II, g = Tro^iH/Lrestj so to 
keep a and E constant across the various-sized rods (so 
that each rod has the same cross-section and is made 
of the same materials), we take the values of g to be 
g = 10.00,8.34,7.15, corresponding to the three choices 
of length. In addition, in accordance with Eqs. (2.1) 
and (2.2), we take k = O.OI/tt^. With these param¬ 
eters, for Lq = 1, we have Tc( 0) = —0.01. We take 
T = t/tc{0) = 1.3,1.5,1.7 and vary T to observe the 
transition. For these three r values, with d = 2 and 
all other elastic parameters corresponding to Lq = Ij 
Tc = 3.78 X 10“®, 6.03 x 10“®, 8.05 x 10“®, satisfying the 
requirement that the persistence length Ip = k/T ~ 10^ 
is much longer than the length of the rod Lq, and, there¬ 
fore, the transverse fluctuations are small. This justifies 
the small tt expansion we make. 


The resulting Ul curves from our MC simulations are 
shown in Fig. 4. Because Tc, as given in Eq. (2.28), de¬ 
pends on the system size Lq through |rc(0)|, it is nec¬ 
essary to shift the Ul curves by theoretical predictions 
of Tc(t) to observe the crossing of the three curves for 
the different system sizes. The crossing of the three Ul 
curves for all three values of r in both two and three di¬ 
mensions verifies our theoretical prediction of the finite- 
temperature buckling transition. 


IV. CONCLUSION AND DISCUSSION 

In this paper, we used both analytic theory and MC 
simulations to investigate the buckling of an extensible 
elastic rod at finite temperature. We find that, in both 
two and three dimensions, buckling is delayed by thermal 
fluctuations, and near the transition, there is a critical 
regime in which the fluctuation correction to the average 
compression force is of order Vt. 

In comparing the two phase diagrams in Figs. I and 2, 
one can observe that the straight-rod phase is more sta¬ 
bilized in three dimensions than in two dimensions. This 
can be intuitively attributed to the fact that in higher di¬ 
mensions, there are an increasing number of transverse, 
soft directions in which segments in the straight rod can 
move compared to the when the rod is buckled. There¬ 
fore, the straight rod is increasingly entropically pro¬ 
tected, as there are a larger number of accessible states. 

Our analytic theory is a perturbative theory that ap¬ 
plies to small fluctuations. This requires that the di¬ 
mensionless temperature T <C 1. This condition can be 
written in terms of the persistence length Ip = k/T as 
Trest/(7r^/p) 1, which is satisfied by stiff {L^est ^ Ip) 
and semiflexible (Trest ~ ^p) polymers. 

At lower temperatures, quantum fluctuations also be¬ 
come important. To make a simple estimate of the tem¬ 
perature scale at which this occurs, we include the kinetic 
energy term 


where p is the linear mass density of the rod. Combining 
this with the potential energy terms in 77, we have a 


.^kinetic — / dXd p 
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T-Tc{L^) T-T,(Lo) T-T,(Lo) 


FIG. 4. Results of Monte Carlo simulations run with f = 1.3,1.5,1.7 in both two and three dimensions. For each value of f 
and in each dimension, three different lengths were simulated. Each data point corresponds to the combined Binder cumulant 
value of ten independent simulations run with identical parameters. The lines connecting the data points explicitly illustrate 
that the Binder cumulant curves do indeed simultaneously cross at, or very close to, the respective critical temperatures. 


phonon energy given by 

hu !^ \/r + (4.2) 

Vp 

Therefore, in addition to thermal fluctuation corrections, 
quantum fluctuations also contribute to the renormaliza¬ 
tion of f and K, moving the critical r to a larger compres¬ 
sion value (in magnitude) even at T = 0. The significance 
of such contributions from quantum fluctuations can be 
estimated by comparing hoj of generic modes with ksT. 
For the simple case of stiff polymers of length 10“® m and 
persistence length 10“^ m, we estimate that the charac¬ 
teristic temperature for huj ~ ksT is T ~ 10“® K, which 
is extremely low. Other systems with higher stiffness or 
shorter lengths may have stronger quantum effects. 

Our result that, in both two and three dimensions, 
the buckling transition is delayed by thermal fluctuations 
contrasts with previous studies of finite-temperature 
buckling of polymers using the inextensible worm-like 
chain model [25, 27]. The extensibility of the rod in 
our model allows for an additional independent quartic- 
order term in the Hamiltonian, and this term plays an 
important role in determining the renormalization of the 
stability of the first mode, leading to the phase diagram 
shown in Fig. 1. 

In addition, extensive recent studies have focused on 
zero-temperature mechanical instability in both ordered 
and disordered systems [32-40], and their behavior at 
finite temperature remain largely unexplored [14, 16- 


19, 41, 42]. Our model provides a clean system which 
exhibit a shifted second-order transition and the results 
can be compared to future studies on finite-temperature 
mechanical instabilities in various systems. 


Appendix A: Deriving the Hamiltonian of the 
extensible rod with fluctuations 


The change in the length of the rod due to thermal 
fluctuations can be expressed in terms of the tt field as 

L-Lo = J dxd ^a/I -I- IttJ^ - 1 ^ . (Al) 

Using this, we can then write the stretching/compression 
elastic energy (2.3) as 


Use — X \-T 

25 Jo 


dxd ( a / 


1 + TT P — 1 


dxd (\/ 


1 -I- TT p — 1 


(A2) 


The bending energy of the rod is given by 


Ub = - / ds ldst(s)l^, 


(A3) 


where s labels the arc length. We assume the bending 
rigidity to be homogeneous along the arc length, given 































































that we are considering the regime where stretching is 
much more energetically costly than bending. Here, t(s) 
is the unit tangent vector at s and |(ist(s)| is the local 
curvature. Ub can also be expressed in terms of Tv{xd)'- 


Ub 



(l + |7r|2)3/2 


(tt • tt')^ 

(l+|7r|2)5/2j ’ 

(A4) 


Appendix B: The T = 0 buckled phase 

As we discussed in the main text, for 4 tc( 0) < r < 
rc(0) at r = 0, the tt = 0 straight state is no longer 
stable. The new stable state has tti ^ 0, and the value 
of TTi is determined by minimizing the total Hamiltonian 
with both 7r2 and tt"* terms. Taking = 0 for all but 
the first mode (|(?| = tt/Lq), the Hamiltonian becomes 


where tt' is shorthand for dn/dxd and we used 

~ _ dr{xd)/dxd _ T^jxd) 

\dr{xd)/dxd\ \T^ixd)\ 

and 


(A5) 


H = Ho + 


T - TcjO) 
4Lo 


Fi 


-|T + 5rc(0) +gLo 


Fi 


(HI) 


The minimum-energy configuration is determined by 


ds = \/l + \TT\‘^dxd- (A6) 

The total Hamiltonian of the rod is a sum of both the 
stretching/compression and the bending contributions, 

H = U,c + Ub. (A7) 

Expanding this Hamiltonian as a series in tt leads to the 
form in Eq. (2.7). 

To obtain the Fourier transform of this Hamiltonian, 
we need to pay special attention to the specific boundary 
conditions of the problem. Here, tt has to be a real-valued 
field, and xj_(xd) (the perpendicular component of r, as 
defined in Eq. 2.5) has to vanish at the two ends, Xd = 0 
and Xd = Lq. This limits the Fourier series of x_L(xd) to 
sm{mrxd/L q) basis functions, and the Fourier series of 
Tv{xd) to cos{mrxd/Lo) basis functions. In order to work 
with the more convenient basis of exponential functions, 
we necessarily extend the rod to Xd G [—Lq, Lq] and limit 
’T^{xd) to be real-valued even functions on this interval 
(correspondingly, x±(xd) is limited to real-valued odd 
functions), so that the value of n^Xd) for —Lq < Xd < 0 
is determined by 

x±(xd) = -X±(-Xd) 



7r(xd) = 7r(-Xd). 

(AS) 

Therefore, 

we can write the Fourier transform as 



2Lo ^ 

(A9) 


7^q= [ dxdT^ixd)e~"'^'"‘^, 

J-Lo 

(AlO) 

with 

Tiir ^ V r ^ 

g=—, nGZ\{0}. 

^0 

(All) 


dH 

Otti 


= 0 , 


(B2) 


where tti denotes the mode corresponding to this 
minimum-energy configuration. Thus, we find that 


Fi| = 



4L2(tc( 0) -r) 


5tc(0) -I- gLo 


(B3) 


Taking the limit of stretching stiffness much greater than 
bending stiffness, gLo ^ |tc( 0)|, |r|, we obtain 


l-^il 


/4Lo(tc( 0) - r) 


(B4) 


This leads to the T = 0 equilibrium buckled configuration 


l*-L(a;d)| 


1- I ^Lo{tc{Q) - t) 
TT V ff 



(B5) 


In two dimensions, where xj_ is simply a number, there 
are two degenerate equilibrium buckled configurations 
corresponding to ±|x_L(xd)|. In three dimensions, how¬ 
ever, there are an infinite number, consistent with a U(l) 
symmetry corresponding to rotation about the x^-axis. 

The energy of this equilibrium buckled configuration 
is 


_ F _ (Te(0) - r)2 _ 2rTe(0) - Tc(0)2 
2 g 2g 2g 


indicating a constant force at T = 0 in the buckled phase 


f = 9 


dH 

Ih 


Tc(0). 


(B7) 


Because '^{xd) is real and even, we have constraints on 
TTq that 


Appendix C: Integrating out fluctuations and 
obtaining the Landau free energy 


TTq = 7r_q = TT*. (AI2) 

Therefore, positive and negative q values do not consti¬ 
tute independent modes. 


In this section, we expand the Hamiltonian in terms 
of TT^ and TT^ and perform the calculation of integrating 
out 7r>. 
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It is clear that 7r< and 7r> are decoupled in the 
quadratic Hamiltonian because they are of different mo¬ 
menta and, therefore, orthogonal, so 


H 2 = H2{7V<) + H2{7T>). 


(Cl) 


In the quartic-order Hamiltonian H 4 , on the other hand, 
they are coupled. 

The partition function of the rod can be written as 

2 ; = J g-y(Ho-l-ff2(7r<)-|-H2(7r>)-|-ff4(7r<.7r>)) 




(C2) 


where Zq is defined in Eq. (2.33) and 




^0 J 


Following a cumulant expansion, we can then write 


,-ff4(7r<,7r>)/T\ 


/; 




Since we are ultimately trying to deduce the effect of 
thermal fluctuations on the stability threshold, we are 
interested in the corrections to the quadratic terms in 
|7r<|. In the straight phase, {'Ka'^b) ~ Tdab, meaning that 
(H4) ~ T|7r<p will provide an 0{T) correction to the 
quadratic-order coefficients, while terms from {Hl)/T 
will result in an 0(T^) correction. Since we are do¬ 
ing a perturbative expansion in small fluctuations, which 
necessitates small temperatures, we need only calculate 


The explicit form of is given in Eq. (2.10), and here we replace tt by tt^ - 1 - tt^. Expanding each term in H 4 out, 
we have 


Lq 4,/2 

J dXd (k|^)> = -^ J dXd {T^aT^a'^bT^b)> 


-Ljl 

Lq 

0 


^^b ^b + i'^b '^b)>+ T^b i^a T^b ) > 


+ '^b ^b (TTa TTa ) > + (’’"a ’’"a ^^b ^b ) > 


^0 r > > 

J dxd |7r<|^-h2(d-l)|7r<|2^^i^^Go5-h4|7r<|2^^i^^Go,-h(|7r>|4); 


0 

Lo 


dxd 


1 


|7r 'I ‘ -I- 2{d + l)|7r^|^ (2Lo)^ ^ 


(C5) 

(C6) 

(C7) 

(C8) 


3k 

'T 


Lq Lq 

J dxd {\n\'^\d^^TT\'^)> =-^ J dxd 


\n<\^\d,,n<\^ + id-l)\7r< 


(2Lo)= 




Oq 


-I- (d — 1)|94; 


.<|2 


(2Lo)- 


J 2 ^ 0 q + (k>n 54 ,^ 7 r>| 2 ); 


(C9) 


Lq Lq 

j dxd ((tt • d^^TTf)> = - ^ J dxd 


(7r< • d^^TV<)‘^ + |7r<|2- 


(2io)= 




Oq 




<|2 


(2Lo)- 


^Go9 + ((7r> 


(CIO) 
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and 


-^0 -^0 

jj dxddx'^ {\Tv{xd)\^\Tv{x'd)\^)> = ^ JJ dxddx'd 


7r<(xd)| \Tv<ix'd)\ +(|7r>(xd)| |7r>(xrf)| )> 


^0 > 

+ |y dxd2{d-l)\n<{xd)fj^^'^Gog. 


(Cll) 


In these equations, 


Go,= 


dLpT 

T + Kq^ 


Using the notation of Eq. (2.23), we can also write 


4L2 


Y,Gog=2fA{f). 

q 


(C12) 


(C13) 


Feynman diagrams corresponding to these terms are included in Fig. 5. 

As mentioned previously, we are interested in extracting the contribution to the coefficients of the quadratic-order 
TT^ terms from Collecting terms, and defining renormalized elastic parameters f and k as the modified 

coefficients, we find 


1 1 
T = T + O + 1)^ + - 1)5^0] 7 dUA 2 


{‘2L^f 


l + -(d-l) 


J (2io)^ 




Og? 


k = K{ \ — 


l + -{d-l) 


(2Ao)‘ 




Og 


(CM) 

(C15) 


> 

The ^ 2 LoY S Q^Goq term appears to have an ultraviolet divergence, but it actually vanishes. This is because it 
9 

originates from quartic-order terms in the bending energy (see Eq. (2.10)) where the spatial derivatives are on the 
legs that combine to form the loops in the Feynman diagrams. This corresponds to a factor of which is the 

leading-order term in the gradient expansion of the difference in orientation between neighboring segments on the rod. 
We can show this by restoring the full form of this factor for a segmented rod, J2xd — '^ixd + , where (q is 

the projected length of each segment, and writing it in momentum space. Doing so, we obtain a factor of 1 — coslgio) 
rather than only the leading-order term q^. Here, q = mr/Nio, so that qio = nir/N. Taking the continuum limit 

> 

where N ^ 00 , the factor l — cos{q£o) is highly oscillatory and the thus the whole expression, S[l~cos((7£o)]G'og 

9 

vanishes. 

Simplifying these equations, we obtain the expressions for the renormalized elastic parameters, f and k, in Eq. (2.22). 


Appendix D: Deriving the effective force and in terms of tti: 


In this section, we derive the effective force as per- 
scribed in Eq. (2.30). Starting from Eq. (2.20) and build¬ 
ing on the calculations of App. C, we have that 

E<(7r<) = i7o + i^2(7r<) + (i/4)>-TlnZo> (Dl) ^nd 

= Ho + F< (7r<) + E4< (7r<) - T In Z> +0{T^), 

with defined as in Eq. (2.21), and con¬ 

taining terms quartic order in tt^. The 0{T^) terms, 
arising from 4-point correlation functions of tt^, can be 
discarded. It is more convenient, going forward, to write 


T + 

F< = —^\7v,f = hT\7v,f (D2) 

4i^0 

< -^T + 5Tc{0)+gLo 4 

=-32Ll- 

«^kil" = &4r|7ri|4, (D3) 
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FIG. 5. Feynman diagrams corresponding to terms in {// 4 )>. The diagrams are systematically divided into rows and columns: 
each row is associated with a single elastic parameter that is the coefficient of the originating term in while the columns 

specify which elastic parameter’s renormalization the diagrams contribute to. Namely, the first column presents the basic vertex 
diagrams, while the second and third columns list those that renormalize r and k, respectively. Each external leg corresponds 
to a TT^ field, and slashes denote spatial derivatives with respect to Xd- The various internal lines differentiate between the 
interactions and are used for index bookkeeping. 


where 


62 = 
64 = 


f k-. 


ATLo 

9 

32TL2 ■ 


(D4) 


Since the first two terms are independent of tti, we can 
easily obtain an expression for the free energy, 

/ OO 

-00 

(D6) 

We now proceed to compute the latter two terms in this 
expression. 

First, 


Once again, we have taken the limit of the stretching 
stiffness much stronger than the bending stiffness in sim¬ 
plifying the expression for . Thus, the Landau free 
energy becomes 


-TlnZo> = -Tin / 


=-Ti„n 


AttTLq 


T + 


4^TLo 


^In 


T -b 


(D7) 


F< = Ho- T\iiZ> + 62F|7ri|2 -b &4T|7ri|4. (D5) 


2 
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Second, we need to evaluate the integral 


^1 




straightforward to see this by plugging tti - as given in 
Eq. (B3) - into H. 

Next, we put the terms together and derive the effec¬ 
tive force. Following Eq. (2.30), 


This integral can be evaluated using the (d — 1)- 
dimensional spherical coordinates; we find 


^1 





r 



iFi 


/ d — 1 1 c\ 
4 ’ 2’ 4/ 


Tv^r 





/d + 1 3 c\ 


(D9) 


= "tT ^^^^5T'|'rc(0)|^(f)-h/i, (D15) 

where fi is from the In Zi term and can be expanded in 
the various limits. 

In the critical regime, we use Eq. (Dll) and find 


where ^d-i is the solid angle subtended by the {d — 1 )- 
dimensional hypersphere, the dimensionless number 



and iFi represents the Kummer confluent hypergeomet¬ 
ric function. The — sign in Eq. (D9) applies to 62 > 0, 
which is the straight phase, whereas the -I- sign corre¬ 
sponds to 62 < 0 , the buckled phase. 

To better understand the expression in Eq. (D9), we 
expand it in different regimes. The behavior of the iFi 
function takes different limits for c <C 1 (close to the 
transition - the critical regime) and c ;g> 1 (far from the 
transition). The boundary between these two regimes, 
determined by c ~ 1 , is indicated by the dashed curves 
in Figs. 1 and 2. For c <C 1, we have 





^d-iT 



1 - f 


^2 


d-l\b 2 ? 

8 bi y 
(Dll) 


X (l-(d-l)gT^'(f)), (D16) 

where A'{f) is the derivative of .4(f) with respect to f. 
Deep in the straight phase, we use Eq. (D12) to obtain 

/i = ^ .2 (1 - . (D17) 

2 T + K(7r/Lo) 

while deep in the buckled phase, Eq. (D13) gives us 

r ( ~ , 3-ci gT \ 

V ^ ^7.2+ 2 f + k{T:/L,Y) 

x{l-{d-l)-gf^{f)). (D18) 

In the latter regimes, where c 1, it turns out that 
gT ^ 1; therefore, we can write the expressions for fi 
to 0{gT). The complete force expressions then simply 
become a leading-order term plus an 0{gT) correction. 
Specihcally, in the straight phase. 



where f = r(^^)/r(^^). For c I, the asymp¬ 
totic expressions depend on the phase of the rod: for 
the straight phase (&2 > 0 ), 




while for the buckled phase (62 < 0 ), 

a/tT <j-3 d-2 

Zi = 1621^64 ^ , 

2 2 




(D12) 


(D13) 


The expressions for Z\ in the c ;g> 1 regime yield simple 
expressions when specializing to d = 2 and d = 3, so it is 
useful to explicitly list them: 


Zi = 


\fK~ib 2 62 > 0 and d = 2 , 

'n jbi 62 > 0 and d = 3, 

^ 27 r/| 52 | 62 < 0 and d = 2 , 

jbd ^2 < 0 and d = 3. 


(DM) 


The factor in the latter two equations comes from 
the finite expectation value of tti when 62 < 0. It is 


z i — r 


so that 


f = r+ ^^-^gT\Tc{0)\ 
and in the buckled phase, 


A{t) -k 


1 — T 


(DI9) 


(D20) 


fi = Tc{0) - T - gT\Tci0)\ 

F {d — I)M(f) F 


(d-I)M'(f)(f-I) 
3-d I 


2 f - 1 


, (D21) 


so that 


/ = 'rc(O) - gT\Tc 


(d-l)M'(f)(f-l) 


d- I ,, , 3-d 1 


(D22) 
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Notice the major difference between the two final expres¬ 
sions for the effective force: deep in the straight phase, 
the force is just the original/unmodified compression 
with a small 0{gT) correction; on the other hand, deep 
in the buckled phase, the force is the zero-temperature 
critical compression with a small correction of the same 
order. 

The critical regime, however, is not constrained to only 
small values of gT, so a similar expansion cannot be 
made everywhere; therefore, we further divide this regime 
into two limiting cases. In the region where gT <C 1, 
Eq. (D16) becomes 

h = + (t - r,(r)) ^ ~ ^ , (D23) 

where we discard all corrections of 0{gT) and also note 
thatr —Tc(0) = r — rc(r)-|-0(gT), using Eq. (2.29). Fur¬ 
thermore, in this regime, r will deviate minimally from 


Tc{T); therefore, we can simply take /i « T^y2gT - which 
is indeed the value of fi on the transition curve - as a rea¬ 
sonable approximation for the entire critical region (for 
gT ^1). Thus, the total force in this regime is 

/ = T + (D24) 

which indicates an 0 {VT) correction to the force in the 
critical regime. 

Finally, when gT ~ 0(1) in the critical regime, the fi 
contribution to the total force is suppressed, as 1 — (d — 
l)gTA'{f) « 0. In this case, 

f = T+‘^^-^gf\Tc{0)\A{f), (D25) 

and there is, once again, an 0{gT) correction to the com¬ 
pression. 
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